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Abstract. Detailed models of a Cen A and B based on new seismological data for a Cen B by Carrier & Bourban 
ll^03ii have been computed using the Geneva evolution code including atomic diffusion. Taking into account the 
numerous observational constraints now available for the a Cen system, we find a stellar model which is in good 
agreement with the astrometric, photometric, spectroscopic and asteroseismic data. The global parameters of 
the a Cen system are now firmly constrained to an age of t = 6.52 ± 0.30 Gyr, an initial helium mass fraction 
Yi = 0.275 ± 0.010 and an initial metallicity {Z/X)i = 0.0434 ± 0.0020. Thanks to these numerous observational 
constraints, we confirm that the mixing-length parameter a of the B component is larger than the one of the A 
component, as already suggested by many authors (Noels et al. 119911 Fernandes & Neuforge ll995l and Guenther & 
Demaraue l2000ll : ae is about 8 % larger than oa (oa = 1.83±0.10 and ae = 1.97±0.10). Moreover, we show that 
asteroseismic measurements enable to determine the radii of both stars with a very high precision (errors smaller 
than 0.3%). The radii deduced from seismological data are compatible with the new interferometric results of 
Kervella et al. 120031 even if they are slightly larger than the interferometric radii (differences smaller than 1 %). 
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1. Introduction 

■ The a Cen system is the ideal target to test our knowledge 
, of stellar physics in solar-like stars, due to its proxim- 
ity which allows precise determination of the fundamental 
, parameters of the two stars. As a result, numerous the- 
oretical analysis of a Cen A and B have already been 
performed. 

Flannery & Ayres H1978|l were the first to compute 
models of the a Cen system. They found that stellar 
models constructed by assuming a solar composition for 
a Cen A and B were not able to reproduce the astromet- 
ric and photometric data. They concluded that the a Cen 
system is more metal-rich than the Sun {ZaCcn/ZQ ^ 2). 

Following the report of a possible detection of p-mode 
oscillations on a Cen A (Fossat et al. I1984f) . Demarque 
et al. I|1986fl investigated the astrometric and physical 
properties of the a Cen system. Firstly, they deduced the 
masses of the components (1.09 ±0.01 and 0.90 ± 0.01 
for a Cen A and B respectively) from orbital measure- 
ments and ground-based parallax. Secondly, they com- 
puted stellar models of a Cen A in order to confront the- 
oretical predicted p-mode frequencies to observed ones. 
They concluded that the observed frequencies were incon- 
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sistent with p-mode spectra constructed from standard 
theoretical models. 

Noels et al. (|1991|l introduced a calibration procedure 
which is based on the simple principle that the four observ- 
ables (two luminosities and two effective temperatures) 
will determine the four unknowns of the modelisation (Z, 
Y , age and mixing-length parameter a, which is assumed 
to be identical for both stars). Using the masses deduced 
by Demarque et al. iP^SSjl . they found Z = 0.04, Y = 0.32, 
a — 1.6 and an age of 5 Gyr. 

Contrary to Noels et al. l(TMT|l . Edmonds et al. l(TM^ 
relaxed the hypothesis of a unique mixing-length param- 
eter for both stars by considering the observed metallic- 
ity as an additional constraint. They were also the first 
to include the effects of helium diffusion. They found 
Y = 0.300 ± 0.005, an age of 4.6 ± 0.4 Gyr, and a mixing- 
length parameter for a Cen A slightly smaller than for 
a Cen B, but did not include an uncertainty analysis to 
determine if this difference was significant. Edmonds et 
al. also computed the theoretical p-mode frequencies of 
a Cen A and B; they predicted a mean large spacing (for 
^ = 0) of 108 /iHz and 179 /iHz and a mean small spacing 
of 6.2 /xHz and 12.6 /iHz for a Cen A and B respectively. 

Neuforge ((1993|l revisited the study by Noels et 
al. H1991|l using OPAL opacities. She obtained a so- 
lution which favours a high value for the metallicity 
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{ZaCcn/ Zq ^ 2) and a unique convection parameter for 
both stars. 

The a Cen system was also studied in order to test 
the modelisation of convection. Lydon et al. ([1993) con- 
structed a series of models of a Cen A and B for the 
purpose of testing the effects of convection modeling both 
by means of the mixing-length theory and by means of 
parametrization of energy fluxes based upon numerical 
simulations of turbulent compressible convection. They 
found that their formulation of convection produced mod- 
els with theoretical radii compatible with the observed 
ones. They were thus able to correctly modelize a star us- 
ing a formulation of convection which does not include an 
adjustable free parameter to determine the radius. 

Fernandes & Neuforge H1995|l also studied the a Cen 
system to test stellar models based upon the mixing- 
length theory of convection and models using the Canuto 
and Mazzitelli formulation (Canuto & Mazzitelli 119911 
I1992II . Their calibration suggested that the mixing-length 
parameter for a Cen B is larger that the one for a Cen A 
if the mass fraction of heavy elements Z is smaller than 
0.038; if Z is larger than this threshold, the two mixing- 
length parameters become very similar. However, because 
of the uncertainties in the observational constraints, these 
differences could not be firmly established. 

Kim H1999|l calibrated the a Cen system using the 
observational constraint [Z / X] of Neuforge- Verheecke & 
Magain H1997|l . His study, which took into account the he- 
lium diffusion, gave an age of about 5.4 Gyr and showed 
that both stars may have the same mixing-length ratio 
(1.6 ~ 1.7). Kim also computed the p-mode frequencies 
of the a Cen system and predicted a mean large spacing 
of 104 and 171 /xHz for a Cen A and B respectively. 

The different calibrations mentioned previously used 
the masses derived by Demarque et al. (1986J: 1.09 ±0.01 
and 0.90 ± 0.01 for a Cen A and B respectively. 
Pourbaix et al. Ijl999 jl performed a simultaneous least- 
squares adjustment of all visual and spectroscopic obser- 
vations of the a Cen system with precise radial velocities 
measurements. They derived new consistent values of the 
orbital parallax, sum of masses, mass ratio and individ- 
ual masses. They found masses of 1.16 ± 0.031 and 
0.97 ± 0.032 Mq for a Cen A and B respectively. The dif- 
ferences between the masses derived by Pourbaix et al. 
H1999|l and the ones obtained by Demarque et al. ( 1986| 
mainly result from the lower parallax found by Pourbaix et 
al. (737.0±2.6 mas instead of 750.6±4.6 mas for Demarque 
et al.). In order to investigate the consequences of these 
new masses, Pourbaix et al. constructed new models of 
the a Cen system and found a significantly smaller age 
than previous estimates (2.7 Gyr). 

Guenther & Demarque H2()()0|l calculated detailed mod- 
els of a Cen A and B based on three different paral- 
laxes: the parallax from Hipparcos (742.12 ± 1.40 mas), 
the parallax of 750.6 ± 4.6 mas obtained by Demarque et 
al. 1,198H| and the value of 747.1 ± 1.2 mas determined 
by Soderhjelm (^99). Using the orbital data from Heintz 
H1982|l and the mass ratio of Kamper & Wesselink H1978|l 



they thus deduced different masses of a Cen A and B 
for each parallax: Ma = 1.1238 ± 0.008 Mq and Mb = 
0.9344 ±0.007Mq for the parallax from Hipparcos, Ma = 
1.0844 ± 0.008 Mq and Mb = 0.9017 ± 0.007 Mq for the 
parallax of Demarque et al., and Ma = 1.1O15±O.OO8M0 
and Mb = 0.9159 ± 0.007 A^q for the parallax deter- 
mined by Soderhjelm. For each pair of masses, Guenther 
& Demarque calculated models of a Cen A and B includ- 
ing helium and heavy-element diffusion which had to re- 
produce the observed effective temperatures, luminosities 
(which are slightly different for each pair of masses due 
to the different parallaxes) and surface metallicities. They 
found that self-consistent models of the a Cen system 
could be produced for each pair of masses and thus that 
the parallax from Demarque et al. as well as the Hipparcos 
parallax could not be ruled out, because observational 
uncertainties in other parameters, such as composition, 
dominated the uncertainties. To investigate in details the 
effect of uncertainties in mass, luminosity, effective tem- 
perature, helium aboundance and metallicity on the mod- 
els, they chose the models based on the most recent de- 
termination of the parallax (the parallax of Soderhjelm). 
For these models, p-mode frequencies were also computed. 
Guenther & Demarque found an initial helium mass frac- 
tion YzAMS — 0.28 and an age of the system which de- 
pends critically on whether or not a Cen A has a convec- 
tive core. If it does, the system has an age of 7.6 ±0.8 Gyr. 
If a Cen A does not have a convective core (correspond- 
ing to ZzAMS < 0.3), a Cen A and B are 6.8 ± 0.8 Gyr 
old. They also found that the mixing-length parameter 
of a Cen A is about 10% smaller that the one of the B 
component. However, they pointed out that the effect of 
composition and surface temperature uncertainties on a 
is greater than this difference. Concerning the pulsation 
analysis of their models, they concluded that the large 
spacing provides a very precise means of determining the 
radius of the star, uncontaminated to any significant de- 
gree by uncertainties in the star's composition. They also 
pointed out that, because a Cen A and B have the same 
age and initial composition, the small spacing may be use- 
ful as an age indicator. They predicted an average large 
and small spacing of 101 ± 3 /iHz and 4.6 ± 0.4 /iHz for 
a Cen A, and 173 ± 6 /iHz and 15 ± 1 /iHz for a Cen B. 

Morel et al. I|2000l calculated detailed evolutionary 
models of the a Cen system based on the new masses 
of 1.16 ± 0.031 Mq and 0.97 ± 0.032 M© for a Cen A and 
B respectively (Pourbaix et al. ll999|l . Contrarily to previ- 
ous works they also chose to use the spectroscopic grav- 
ities instead of the luminosity derived from the photom- 
etry, bolometric correction and parallax. The evolution- 
ary code they used included microscopic diffusion. With 
the mixing-length theory of convection. Morel et al. found 
that the a Cen system has an age of 2.71 Gyr (in very good 
agreement with Pourbaix et al. I1999|l . an initial chemical 
composition Yi = 0.284 and {Z/X)i = 0.0443, and val- 
ues of the convection parameters which are almost the 
same for both stars. Including overshooting of convective 
cores, they found a larger age of the system (3.53 Gyr), 
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an initial chemical composition almost the same as with- 
out overshooting {Y, = 0.279 and {Z/X)-, = 0.0450), 
and mixing-length parameters almost identical for both 
components. Morel et al. also performed a calibration of 
the a Cen system using the Canuto and Mazzitelli con- 
vection theory. They found that this formulation of the 
convection changed the results of the calibration. Indeed, 
the age of the a Cen system is larger with the Canuto 
and Mazzitelli convection theory than with the mixing- 
length theory (4.086 Gyr instead of 2.71 Gyr). The initial 
composition derived using the Canuto and Mazzitelli for- 
mulation is Yi = 0.271 and {Z/X)-, = 0.0450. The con- 
vection parameters of both components were also almost 
identical and close to unity. Morel et al. also calculated 
models using the observational constraints adopted by 
Guenther & Demarque H20[)0|l . This calibration gave a 
smaller age (5.64 Gyr) and a larger initial helium mass 
fraction (li = 0.300) than obtained by Guenther & 
Demarque l)2000|l . Finally, Morel et al. calculated the p- 
mode frequencies of their models based on the masses of 
Pourbaix et al. H1999(l . They predicted large and small 
spacings of about 107 — 108 /iHz and 7.5 — 9.1 /iHz for 
a Cen A and 154 — 157/iHz and 12 /iHz for a Cen B. 

The detection and identification of p-modes in 
a Cen A by Bouchy & Carrier (j2()02|l brought strong ad- 
ditional constraints. These measurements were not fully 
in agreement with the theoretical calibrations quoted 
above, which were based on non-asteroseismic observ- 
ables. Moreover, Pourbaix et al. (2002J studied in detail 
the masses of the a Cen system. Thanks to the accu- 
rate estimate of the parallax fSoderhielm I1999jl and new 
radial velocities, they determined very precise masses of 
1.105±0.0070Mo and 0.934±0.0061Mq for a Cen A and 
B respectively. 

Consequently, new calibrations of the a Cen system 
were performed by Thevenin et al. I|2002|l and Thoul et 
al. H2003|l in order to find a model reproducing the as- 
teroseismic observations of the A component as well as 
the new masses of Pourbaix et al. H2002|l . The analysis 
of these two theoretical groups led to different results. 
Thevenin et al. found a model with an age of 4.85 Gyr 
which is not able to reproduce the astrometric mass of 
Pourbaix et al. H20()2|l of a Cen B, whereas Thoul et al. 
proposed an older model (6.41 Gyr) which matches astro- 
metric masses of both stars. The main difference between 
these two analysis results from the inclusion of atomic dif- 
fusion. Indeed, the evolution code used by Thoul et al. 
does not include the diffusion of helium and other heavy 
elements which cannot be neglected in order to obtain ac- 
curate stellar models of stars in the mass range of a Cen A 
and B. On the other hand, Thevenin et al. included the 
atomic diffusion in their evolution code and assumed the 
mixing-length parameter of the two components to be 
the same. Moreover, Thevenin et al. used the luminos- 
ity derived from the photometry, bolometric correction 
and parallax contrary to Thoul et al. who used the spec- 
troscopic surface gravities to constrain the system. Both 
studies mentioned that, even if the p-mode frequencies 



of a Cen A bring strong constraints for the calibration of 
the a Cen system, asteroseismic measurements of a Cen B 
are needed to unambiguously determine the fundamental 
parameters of a Cen A and B. 

Recently, Carrier & Bourban 1)2003(1 detected solar- 
like oscillations in a Cen B with the Coralie echelle 
spectrograph and identified twelve individual frequencies. 
Thanks to these new asteroseismic measurements, there 
are now enough observational constraints to determine an 
accurate model for the a Cen system. Moreover, Kervella 
et al. I;2003,l recently measured the angular diameters of 
a Cen A and B. These results enable us to test the con- 
sistency between interferometric and asteroseismic obser- 
vations. 

The whole observational constraints available for the 
a Cen system are presented in Sect. 2. The physics of the 
stellar models and the calibration method are described 
in Sect. 3, while the results are given in Sect. 4. Finally 
Sect. 5 contains the comparison of our results with previ- 
ous studies and the conclusion is given in Sect. 6. 

2. Observational constraints 

2.1. Astrometric data 

The parallax and the orbital elements of the a Cen bi- 
nary system are difhcult to determine accurately due to 
the high luminosity and large separation of the stellar 
components. This results in various values for the par- 
allax published in the literature (see the introduction and 
Table 10 of Kervella et al. In this work we con- 

sider the accurate parallax obtained by Soderhjelm l|1999|l 
that combines Hipparcos and ground-based observations, 
747.1 ±1.2 mas. 

Pourbaix et al. H2002fl improved the precision of the 
orbital parameters of the a Cen system by adding new 
radial velocity measurements of a Cen A and B ob- 
tained in the framework of the Anglo-Australian Planet 
Search programme as well as in the CORALIE programme 
to those by Endl et al. H2001|l . Adopting the parallax 
of Soderhjelm H1999|l . they determined precise masses of 
1.105 ± 0.0070 Mq and 0.934 ± 0.0061 Mq for a Cen A 
and B respectively. Note that the mass ratio they deter- 
mined is compatible with the one obtained by Kamper & 
Wesselink H1978|l . These new masses constitute the most 
recent and accurate values now available for a Cen A and 
B. Consequently, we consider them as true observables 
which have to be reproduced by a consistent model of the 
a Cen system. 

2.2. Effective temperatures and chemical composition 

Many spectroscopic measurements of both components of 
the system can be found in the literature. A summary of 
published values is given in Table 3 of Morel et al. H2000|l . 

Concerning the effective temperatures, Thevenin et al. 
(2002) and Thoul et al. plH^ used the same value for 
a Cen B. In this work, we considered the same effective 
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Table 1. Geneva photometric measurements of a Cen A 
and B. P corresponds to the weight which varies from 
to 4 according to the quahty of the nights (4 for the best 
nights and for the worst). 



HJD - 2 430 000 Vk 


Vb 


P 


946.718 




1.296 


1 


946.720 


0.006 


1.361 


1 


1056.543 


-0.023 


1.357 


1 


1063.558 


0.003 


1.331 


1 


1070.518 


0.006 


1.344 


1 


1324.808 


0.001 




1 


1688.823 


-0.005 


1.324 


3 



temperature of a Cen B of 5260 ± 50 K. For the A com- 
ponent, we note a small discrepancy between the value of 
Neuforge-Verheecke & Magain H1997|l used by Thoul et al. 
and the temperature derived by Morel et al. lj^OO|) and 
used by Thevenin et al. We thus adopted an effective tem- 
perature of 58 10 ± 50 K for a Cen A in order to encompass 
the intervals of temperatures used by Thevenin et al. and 
Thoul et al. 

For the metallicities, we adopted [Fc/HJa = 0.22 ± 
0.05 dex for a Cen A and [Fe/H]B = 0.24 ± 0.05 dex 
for a Cen B. These values lie between the ones used by 
Thevenin et al. and those adopted by Thoul et al, with 
larger error boxes which appear to us more realistic in 
view of the different values found in the literature. 

2.3. Luminosities 

From 1978 to 1981, both components of the a Cen sys- 
tem have been measured in the Geneva photometric sys- 
tem fGolav [l980|l with the photoelectric photometer P7 
(Burnet & Rufener ri979|l installed on the 40 cm Swiss tele- 
scope in La Silla (ESQ, Chile). Six measurements were 
obtained for both stars (see Table ^ using a mask on the 
telescope in order to reduce the flux of these stars; the 
photometric reduction procedure is described by Rufener 
(1964, 1985). Combining the mean magnitudes (Va) = 
-0.003 ± 0.006 and {Vb) = 1.333 ± 0.014 mag, the par- 
allax of Soderhjelm H1999|l . the solar absolute bolomet- 
ric magnitude Mboi, © = 4.746 (Lejeune et al. 1998) and 
the bolometric corrections from Flower's (1996) BCp^ = 
-0.074 ± 0.009 and BC'b = -0.207 ± 0.017 mag deter- 
mined from the effective temperatures, luminosities are 
calculated and have values of La = 1.522 ± 0.030^0 and 
Lb = 0.503 ± 0.020 Lq. 

2.4. Angular diameters 

Recently Kervella et al. ()2()()3(l measured the angular di- 
ameters of a Cen A and B using the VINCI instrument 
installed at ESO's VLT Interferometer. They found limb 
darkened angular diameters 9a — 8.511 ± 0.020 mas and 
9b = 6.001 ± 0.034 mas for the A and B component re- 
spectively. The radius of each star can then be deduced by 



Table 2. Observational constraints for a Cen A and B. 
References: (1) Soderhjelm (1999), (2) Pourbaix et al. 
(|2002(l . (3) this paper, (4) derived from the other observa- 
tional measurements (see text), (5) Kervella et al. H2003|l . 
(6) Bouchy & Carrier (^2002) and (7) Carrier & Bourban 

(12003). 





a Cen A 


a Cen B 


References 


n [mas] 


747.1 ± 1.2 


(1) 


M/Mq 


1.105 ±0.0070 


0.934 ± 0.0061 


(2) 


V [mag] 


-0.003 ±0.006 


1.333 ±0.014 


(3) 


L/Le 


1.522 ±0.030 


0.503 ± 0.020 


(4) 


Teff [K] 


5810 ± 50 


5260 ± 50 


(3) 


[Fe/H], 


0.22 ±0.05 


0.24 ±0.05 


(3) 


9 [mas] 


8.511 ±0.020 


6.001 ± 0.034 


(5) 


R/Rq 


1.224 ±0.003 


0.863 ± 0.005 


(5) 




105.5 ±0.1 


161.1 ±0.1 


(6), (7) 


Suo2 [mHz] 


5.6 ±0.7 


8.7 ±0.8 


(6), (7) 



using the following simple relation between the angular di- 
ameter 9 (in mas), the radius of the star R (in Rq) and the 
parallax tt (in mas) : 9 = 9.305T0~^(2i?)7r. By taking into 
account the error on the parallax of Soderhjelm H1999II . 
we obtained a radius of 1.224 ± 0.003 i?© for a Cen A and 
0.863 ± 0.005 i?0 for a Cen B. 

2.5. Asteroseismic constraints 

Solar-like oscillations in a Cen A have been detected by 
Bouchy & Carrier H2002|l with the Coralie echelle spec- 
trograph. Twenty-eight oscillation frequencies were iden- 
tified in the power spectrum between 1.8 and 2.9 mHz 
with amplitudes in the range 12 to 44 cm s~^. The ra- 
dial orders of these modes with angular degrees ? < 2 lie 
between 15 and 25. According to the asymptotic theory 
('Tassoul ll980f) the oscillation spectrum can be character- 
ized by two frequency separations, the large and the small 
spacing. The large frequency spacing Ai/i = — Vn-i,i 
corresponds to differences between frequencies of modes 
with the same angular degree I and consecutive radial 
order. For high frequency oscillation modes, this spacing 
remains approximately constant with a mean value Avq 
proportional to the square root of the star's mean den- 
sity. Thus Avq puts strong constraints on the value of the 
radius of the star and enables us to test the consistency 
between interferometric measurements and asteroseismic 
data. The small spacing 5^1^+2 = Vn+i,i — Vn,i+2 is the dif- 
ference between the frequencies of modes with an angular 
degree I of same parity and with consecutive radial order. 
This small spacing is very sensitive to the structure of the 
core and mainly to its hydrogen content, i.e. to its evolu- 
tionary status. Unlike the large spacing, (5z^;,;+2 can vary 
sensitively with the frequency. This fact has to be taken 
into account when one compares theoretical value of this 
mean small spacing with asteroseismic results. Bouchy & 
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Carrier determined the averaged large and small spacing 
of a Cen A by a least square fit of the asymptotic relation 
with their twenty-eight oscillations frequencies identified. 
They found a mean large spacing Ai^,^ = 105.5 ±0.1 /j,Hz 
and a mean small spacing between I = 2 and I — modes 
Su^2 ~ 5.6±0.7/iHz. Note that their fit also gives the value 
of a third parameter, e, which is sensitive to the reflexion 
conditions near the surface of the star. Given that the ex- 
act values of the frequencies depend on the details of the 
star's atmosphere, where the pulsation is non-adiabatic, 
we will not use this parameter to constrain our stellar 
models. Indeed, a linear shift of a few /iHz between the- 
oretical and observational frequencies is perfectly accept- 
able. 

Recently, Carrier & Bourban (|2()()H|I detected solar- 
like oscillations in a Cen B with the Coralie echelle 
spectrograph. They identified twelve modes between 3 and 
4.6 mHz with amplitudes in the range 8.7 to 13.7cm s^^. 
From these twelve frequencies they deduced a large spac- 
ing Ai'^ — 161.1 ± 0.1 /iHz and a mean small spacing 
Sug = 8.7±0.8/mz. 

All the observational constraints are listed in Table |21 

3. Stellar models 

3.1. Input Physics 

The stellar evolution code used for these computations is 
the Geneva code, described several times in the literature 
(see Meynet & Maeder 2000lfor more details) . We used the 
OPAL opacities, the MHD equation of state (Dappen et al. 
\WM Hummer & Mihalas.1988. Mihalas et al. 1988), the 
NACRE nuclear reaction rates (Angulo et al. 1999) and 
the standard mixing-length formalism for convection. 

Our models have been computed including atomic dif- 
fusion on He, C, N, O, Ne and Mg using the routines 
developed for the Geneva-Toulouse version of our code 
(see for example Richard et al. ll996|l recently updated by 
O. Richard (private communication). The diffusion coef- 
ficients are computed with the prescription by Paquette 
et al. H1986p . We included the diffusion due to the con- 
centration and thermal gradients, but the radiative accel- 
eration was neglected as it is negligible for the structure 
of the low-mass stellar models with extended convective 
envelopes (Turcotte et al. I1998|l . 

3.2. Calibration method 

Basically, the calibration of a binary system consists in 
finding the set of stellar modeling parameters which best 
reproduces all observational data available for both stars. 
For a given stellar mass the characteristics of a stellar 
model (luminosity, effective temperature, surface metal- 
licity, frequencies of oscillation modes, etc.) depend on 
four modeling parameters: the age of the star (noted t 
in the following), the mixing-length parameter a = I /Hp 
for convection and two parameters describing the ini- 
tial chemical composition of the star. For these two pa- 



rameters, we chose the initial helium abundance Yi and 
the initial ratio between the mass fraction of heavy ele- 
ments and hydrogen {Z/X)i. Assuming that this ratio is 
proportional to the abundance ratio [Fe/H], we can di- 
rectly relate (Z/X) to [Fe/H] by using the solar value 
iZ/X)Q = 0.0230 given by Grevesse & Sauval ijTOMjl . 
Thus, any characteristic ^ of a given stellar model has 
the following formal dependences with respect to model- 
ing parameters : A — A{t, a, Y^, {Z/X)-/). 
The binary nature of the system provides three con- 
straints: iA = iB, = and {Z/X)f- = {Z/X)f. 
Consequently, we have to determine a set of five modeling 
parameters (i, oa, Q;b, and {Z/X)\) instead of eight 
(four for each star). Note that we do not assume that the 
mixing-length value is identical for a Cen A and B. This 
assumption was often used in the past (Noels et al. I1991|l 
as an additional constraint needed to close the system, 
because the only observed values were the effective tem- 
peratures, the masses and the luminosities of both stars. 
We no longer need this assumption because we now have 
enough observational data to strongly constrain the mod- 
els. Moreover, one of the applications of asteroseismology 
is to determine whether or not a single mixing-length pa- 
rameter is applicable to all stars regardless of mass, com- 
position and age. 

3.2.1. Non-asteroseismic constraints 

Once assuming that the masses of a Cen A and B are 
those determined by Pourbaix et al. H20()2|l . the determi- 
nation of the set of modeling parameters (t, oa, Q^b, Yi 
and {Z/X)i) leading to the best agreement with the ob- 
servational constraints is made in two steps. First, we only 
consider non-asteroseismic observations and construct a 
grid of models with positions of the two components in the 
HR diagram in agreement with the observational values of 
the luminosities, effective temperatures and radii listed in 
Table 12 The boxes in the HR diagram for a Cen A and B 
are shown in Fig. ^ To construct the grid of models with 
positions of the two components lying in the observational 
boxes, we proceed in the following way: for a given chem- 
ical composition (i.e a given set Yi, {Z/X)i) the mixing- 
length coefficient of each star is adjusted in order to match 
the observational position in the HR diagram. Note that 
because of the uncertainties on the effective temperature, 
different values of the mixing-length parameter a enable 
to match the observational box in the HR diagram. Thus, 
for each component and for a given initial chemical com- 
position, we obtain different models corresponding to dif- 
ferent value of a between the smallest value of a needed 
to match the observational box (denoted amin) and the 
largest value (amax)- The uncertainty of 50 K on the effec- 
tive temperature of both stars leads to a typical difference 
between amin and amax of about 0.1. 

Once the positions in the HR diagram of the two com- 
ponents agree with the observed values, the surface metal- 
licity of the two stars is compared to the observed one. 
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If the surface metallicity of one of the component is out 
of the metaUicity intervals hsted in Table |2 the models 
are rejected and the procedure is repeated with another 
choice of and {Z/ X)^. Note that the surface metallicities 
[Fe/H]s are almost identical for the models with the same 
initial composition and different mixing-length parame- 
ters a. Moreover, the [Fe/H]s of the models are mainly 
sensitive to {Z/X){ and less to Yi. As a result, the val- 
ues of {Z/X)i are directly constrained by the observed 
surface metallicities; we found that the models matching 
the observed metallicities have {Z/X)i ranging from about 
0.038 to 0.048 (i.e. an initial metallicity [Fe/H]i between 
0.22 and 0.32). 

If the metallicities of both stars agree with the ob- 
served values, we then compare their ages. All the mod- 
els of a Cen A with different values of the mixing-length 
parameters between amin and oimax are considered. For 
all these models, whose position in the HR diagram as 
well as the surface metallicity are in agreement with the 
observational values, the smallest (t^j^) and the highest 
ages (imax) are determined. Note that, given the shape of 
the observational box in the HR diagram (see Fig. QJ, 
the youngest model has a = amin and the oldest has 
a = Qfmax- In the same way, t^^^ and t^^^^ are determined 
for the B component. If t^-^^^ > t^^^ or t^^^ < t^^^ the 
models of a Cen A are not compatible in age with the 
models of the B component; they are rejected and the pro- 
cedure is repeated with another choice of Yi and {Z/X)i. 
Otherwise, all the models of a Cen A and B with the same 
age (with a difference smaller than 0.01 Gyr) are consid- 
ered as models of the a Cen system which reproduce all 
the non-asteroseismic constraints. The whole procedure is 
then repeated with a new choice of Yi and {Z/X)i. 

In this way, we obtained a grid of models with various 
sets of modeling parameters {t, aA, ae, Yi and {Z/X)i) 
which satisfied all the non-asteroseismic observational 
constraints, namely the effective temperatures, the lumi- 
nosities, the radii and the surface metallicities. The second 
step in determining the best model of the binary system 
a Cen is to consider the asteroseismic measurements for 
both stars. Note that if no asteroseismic observations were 
available for the a Cen system, we would not be able to 
discriminate between these stellar models. Indeed, aster- 
oseismic observations are needed to differentiate models 
with different internal structures located in the same re- 
gion of the HR diagram, which is absolutely necessary in 
order to determine the age of the system for instance. 

3.2.2. Pulsation analysis 

For each stellar model of the grid constructed as explained 
above, \ow-l p-mode frequencies were calculated using 
the Aarhus adiabatic pulsation package (Christensen- 
Dalseaard 11997(1 . Following the observations, frequencies 
of modes of degree I < 2 were calculated for a radial or- 
der n ranging from 15 to 25 for a Cen A and from 17 
to 27 for a Cen B. The determination of the large and 
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Fig. 1. Observational constraints in the HR diagram for 
a Cen A and B. The dashed lines indicate the boxes de- 
limited by the observed luminosities and effective temper- 
atures (with their respective 1-sigma errors). The con- 
tinuous lines denote the boxes delimited by the observed 
radii with their 1-sigma errors, while the dotted line cor- 
respond to the same observed radii with errors of 2 sigma. 
The positions in the HR diagram of a Cen A and B for 
the Ml model is indicated by the two triangles, while the 
dots correspond to the positions of the M2 model. Note 
that for a Cen B the dot and the triangle are partly su- 
perimposed, because the Ml model is only slightly hotter 
and less luminous than the M2 model. 

small spacings was carried out exactly as in Bouchy & 
Carrier i^Ud^ and Carrier & Bourban For the 

models of a Cen A, we selected the theoretical modes cor- 
responding to the twenty-eight frequencies identified by 
Bouchy & Carrier and fitted the asymptotic relation to 
obtain the mean values Au^ and Si'q2- In the same way, 
we selected the theoretical modes of the a Cen B mod- 
els corresponding to the twelve frequencies identified by 
Carrier & Bourban and fitted the asymptotic relation to 
obtain Avq and 6l'q2- 

Once the asteroseismic characteristics of all the rele- 
vant models were determined, we performed a mini- 
mization in order to deduce the set of parameters {t, qa, 
as, Yi, {Z/X)i) leading to the best agreement with the 
observations. Therefore we defined the Xtot functional 



Xtot 



EE 



A,B 



a 



tlico 



(1) 



where the vectors C contains all the observables for one 
star: 

C = (i/LQ,Teff,i?/i?0,[Fe/H]„Ai^o,'^i^O2) • 

The vector C*''°° contains the theoretical values of these 
observables for the model to be tested, while the values of 
C"*^^ are those listed in Table [3 The vector aC contains 
the errors on these observations which are also given in 
Tabled 

To better test the agreement between models and as- 
teroseismic observations, we defined a second functional 
Xastoro which dircctly compares individual theoretical fre- 
quencies to the observed ones instead of using the mean 
large and small spacings. As mentioned in Sect. 12.51 a 
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linear shift of a few /iHz between theoretical and observa- 
tional frequencies is perfectly acceptable, because the ex- 
act values of the frequencies depend on the details of the 
star's atmosphere, where the pulsation is non-adiabatic. 
To take this fact into account, we defined the mean value 
of the differences between the theoretical and observed 
frequencies : 

1 ^ 

i=l 

where N is the number of observed frequencies (A^a — 28 
for a Cen A and TVb = 12 for a Cen B). The functional 
Xastero ^an then be defined as 




where (ta = 0.46 /iHz and tiB — 0.47 /iHz are the errors on 
the observed frequencies for a Cen A and B respectively. 

The determination of the best set of parameters (i, 
aA, ftB, Yh I was based on the minimization of the 
functional defined in equation which includes four non- 
asteroseismic and two asteroseismic observational con- 
straints for each star. Once the model with the smallest 
Xtot "^as determined, we refined the grid in the vicinity of 
this preliminary solution in order to find the best solution 
which minimizes at the same time x^^^. and Xastero- 

4. Results 

4.1. Model Ml: 1-sigma error on the observed radii 

We first computed a grid of models reproducing all non- 
asteroseismic constraints within their one-sigma error 
boxes as described in Sect. 13.2.11 The positions in the 
HR diagram for a Cen A and B are given by their lu- 
minosities, effective temperatures and radii (see Fig. ^J. 
Of course, these independent measurements are not fully 
consistent with each other; for instance, the radius deter- 
mined by interferometry has not exactly the same value 
as the one deduced from the luminosity and the effective 
temperature of the star. 

Once this grid of models was computed, we performed 
the minimization described above to find the set of 
modeling parameters which best reproduced all obser- 
vational constraints. In this way, we found the solution 
t = 6.50 ± 0.20 Gyr, aA = 1.83 ± 0.10, ae = 1.99 ± 0.10, 
Yi = 0.275 ± 0.010 and {Z/X)i = 0.0435 ± 0.0020. The 
position in the HR diagram of this model of a Cen A and 
B (denoted model Ml in the following) is given in Fig. ^ 
The characteristics of this model are reported in Table |21 
The confidence limits of each modeling parameter given in 
Tableware estimated as the maximum/minimum values 



which fit the observational constraints when the other cal- 
ibration parameters are fixed to their medium value. The 
asteroseismic features of this solution are given in Fig. [5] 
and Fig. O These two figures show that this solution is 
not in complete agreement with the asteroseismic obser- 
vations. Indeed, the theoretical mean large spacing Avq of 
a Cen A and B are respectively 4cr and 6a larger than the 
observed ones. This can be seen in Fig. [21 which exhibits 
the large spacing versus frequency for both stars, but is 
more clearly shown in Fig. |31 where differences between 
calculated and observed frequencies are plotted. In fact, 
the differences between calculated and observed frequen- 
cies increase with frequency. This trend, which is more 
pronounced for a Cen B than for a Cen A, results di- 
rectly from the too large values of the theoretical large 
spacings. The solution which minimized the functional 
implies that the two stars lie in the upper border of their 
respective boxes in radius (see Fig.^. Given that the large 
spacing Av of a star is proportional to the square root of 
its mean density, we find that the radii observed by in- 
terferometry are slightly smaller that the ones needed to 
reproduce the asteroseismic observations of a Cen A and 
B. 
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Fig. 2. Large spacing versus frequency for a Cen A and 
B for the model with 1-sigma error boxes on the observed 
radii (model Ml). The dots indicate the observed values 
of the large spacing with an uncertainty on individual fre- 
quencies estimated to half the time resolution (0.46 and 
0.47 /iHz for a Cen A and B respectively). 



4.2. Model M2: 2-sigma error on the observed radii 

In order to determine a model of a Cen A and B which 
is in better agreement with the asteroseismic measure- 
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Table 3. Models for a Cen A and B. The upper part of the table gives the observational constraints used for the 
calibration. The middle part of the table presents the modeling parameters with their confidence limits, while the 
bottom part presents the global parameters of the two stars. The mean large and small spacings and their respective 
errors are calculated exactly as in Bouchy & Carrier l)2(K)2|l and Carrier & Bourban H2()()H|I . 





Model Ml 


Model M2 




a Cen A 


a Cen B 


a Cen A 


a Cen B 


M/Mq 


1.105 


0.934 


1.105 


0.934 


L/Lq 


1.522 ± 0.030 


0.503 ± 0.020 


1.522 ±0.030 


0.503 ± 0.020 


%B [K] 


5810 ± 50 


5260 ± 50 


5810 ± 50 


5260 ± 50 


R/ Rq 


1.224 ±0.003 


0.863 ± 0.005 


1.224 ±0.006 


0.863 ±0.010 


[Fe/H]e 


0.22 ± 0.05 


0.24 ±0.05 


0.22 ± 0.05 


0.24 ± 0.05 




105.5 ±0.1 


161.1 ±0.1 


105.5 ±0.1 


161.1 ±0.1 




5.6 ±0.7 


8.7 ±0.8 


5.6 ±0.7 


8.7 ±0.8 


t [Gyr] 


6.50 ± 0.20 


6.52 ±0.30 


Q 


1.83 ±0.10 


1.99 ±0.10 


1.83 ±0.10 


1.97 ±0.10 


Yi 


0.275 ±0.010 


0.275 ±0.010 




0.0435 ± 0.0020 


0.0434 ± 0.0020 


L/Lq 


1.497 


0.522 


1.503 


0.523 


Teff [K] 


5769 


5270 


5769 


5266 


R/Rq 


1.227 


0.868 


1.229 


0.870 




0.231 


0.247 


0.231 


0.247 




0.0386 


0.0402 


0.0385 


0.0402 


[Fe/H]e 


0.22 


0.24 


0.22 


0.24 




105.9 ±0.1 


161.7 ±0.1 


105.5 ±0.1 


161.1 ± 0.1 


6^02 [mHz] 


4.6 ±0.6 


10.3 ±0.9 


4.6 ±0.6 


10.2 ±0.8 
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Fig. 3. Differences between calculated and observed fre- 
quencies for the model with 1-sigma error boxes on the ob- 
served radii (model Ml). The systematic shifts for a Cen A 
and B are (£',/}a = — 11.5/iHz and {D^)b — IS.TyitHz (see 
text for more details). 

ments, we decided to repeat the above calibration consid- 
ering larger values for the errors on the observed radii. 
We thus took an error on the observed radii of two sigma 
instead of the one-sigma error considered previously: R — 



1.224 ± 0.006 Rq for a Cen A and i? = 0.863 ± 0.010 i?© 
for a Cen B. For the other observational constraints, we 
kept the same values as the previous ones. These new in- 
tervals in radius defined a larger region in the HR dia- 
gram which is shown in Fig.Q] With these new constraints 
on the radii, we found the solution t = 6.52 ± 0.30 Gyr, 
aA = 1.83 ± 0.10, ae = 1.97 ± 0.10, Yi = 0.275 ± 0.010 
and iZ/X)i = 0.0434 ± 0.0020. The position of a Cen A 
and B in the HR diagram for this model (noted model M2 
in the following) is denoted by a dot in Fig. ^ The char- 
acteristics of this model are given in Table 13 Note that 
the radii of a Cen A and B for this solution are only 1.7 
and 1.4 sigma larger than the observed values. Table|3|and 
Fig-Elclearly show that the mean large spacings of the M2 
model are in perfect accordance with the observed values. 
Indeed, the trend visible for the Ml model in Fig. |31 is 
no longer present in the M2 model. The variations of the 
large and small spacing with the frequency for a Cen A 
and B are given in Fig.[Sland|Hl One can see that the M2 
model perfectly reproduces the observed large spacings. It 
is interesting to notice that asteroseismic observations en- 
able an accurate determination of the radii of both stars. 
Indeed, one can see that the Ml model with a radius of 
1.227 i?0 for a Cen A and 0.868 i?0 for a Cen B is not 
compatible with the observed large spacings, while the M2 
model with slightly larger radii (1.229 i?© for a Cen A and 
0.870 Rq for a Cen B) is in perfect agreement with the ob- 
served large spacings. Consequently, we deduce that the 
asteroseismic measurements enable to determine the radii 
of both stars with a very high precision (errors smaller 
than 0.3%). 
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Concerning the small spacings Si'02, we note that the 
agreement with the observed values is not as good as for 
the large spacings. For a Cen A, the theoretical mean 
small spacing is compatible with the observational mea- 
surement but remains slightly smaller than the observed 
value of 5.6 /iHz. For a Cen B, this is exactly the opposite: 
the theoretical mean small spacing is higher than the ob- 
served one. Since the small spacing Si'q2 is very sensitive to 
the hydrogen content of the core of a main-sequence star, 
i.e. mainly to its age and its initial chemical composition, 
we conclude that the age of about 6.5 Gyr given by the 
calibration constitutes the best compromise to reproduce 
at the same time the small spacing of a Cen A and that 
of a Cen B. Note that the observed mean small spacing 
of a Cen B is not strongly constrained since it is deduced 
from only two points (see Fig. inj. The point at low fre- 
quency is perfectly reproduced by the model, contrary to 
the value at higher frequency which is smaller than pre- 
dicted by the model. As a result, the theoretical decrease 
of i5z/o2 with frequency for a Cen B is smaller than ob- 
served. If this trend for a Cen B can be thought to be 
artificial given that it is deduced from only two observa- 
tional points, it is interesting to mention that this trend 
is also present for a Cen A (see Fig. Isj). Note that for 
a Cen A this discrepancy is always present even for stel- 
lar models with quite different ages (see Fig. 2 of Thevenin 
et al. nna^ and Fig. 4 of Thoul et al. IMT^ and remains 
therefore difficult to explain. 
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Fig. 4. Differences between calculated and observed fre- 
quencies for the model with 2-sigma error boxes on the 
observed radii (model M2). The systematic shifts between 
theoretical and observed frequencies for a Cen A and B 
are {Dy)p^ = — 19.0/iHz and {D^)-^ — 0.7 /iHz (see text 
for more details). 
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Fig. 5. Large and small spacings versus frequency for 
a Cen A for the model with 2-sigma error boxes on the 
observed radii (model M2). The dots indicate the observed 
values of the large and small spacings with an uncertainty 
on individual frequencies estimated to half the time reso- 
lution (0.46 and 0.47 /xHz for a Cen A and B respectively). 
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Fig. 6. Large and small spacings versus frequency for 
a Cen B for the M2 model. The dots indicate the observed 
values of the large and small spacings with an uncertainty 
on individual frequencies estimated to 0.46 and 0.47 /iHz 
for a Cen A and B respectively. 

Previous analysis of the a Cen system made without 
the asteroseismic constraints on the B component dis- 
agreed on the presence of a convective core in a Cen A. 
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Thevenin et al. l|2l)U2fl found a model with a convectivc 
core, while Thoul et al. P003) proposed a model without 
a convective core. Our model of a Cen A does not have 
a convective core. However, it is important to note that 
a Cen A lies very close to the boundary between models 
with and without convective core. Thus small changes in 
the observational constraints adopted to calibrate the sys- 
tem can lead to models of a Cen A with or without a con- 
vective core. Moreover, the asteroseismic measurements of 
a Cen A are not precise enough to allow a direct discrim- 
ination on this characteristic. As can indeed be seen on 
Fig. [3 the spacing dv^ = Vn+i,i=a + Vn,i=a - 2vn,i=i for 
our best model of a Cen A is very similar to the spac- 
ing dvi^i of a model with a convective core (see Fig. 3 of 
Thevenin et al. I2()()2|l . Therefore it cannot be used to dis- 
criminate models with and without a convective core. 
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Fig. 7. 5v{)i spacing versus frequency for a Cen A for the 
M2 model. The uncertainty on the observed individual 
frequencies is estimated to 0.46 /iHz. 



Finally, we compare the theoretical p-mode frequen- 
cies of our best model (model M2) to the observed fre- 
quencies by plotting the echelle diagrams of the two stars. 
Fig. |S1 shows the echelle diagram of a Cen A, while that 
of a Cen B is given in Fig. O In these two figures, the 
systematic differences between theoretical and observed 
frequencies {{D^)x = —19.0 //Hz and {D^)b = 0.7 fiRz 
for a Cen A and B) have been taken into account. These 
two figures show the good agreement between our model 
of the a Cen system and the asteroseismic observations. 
The theoretical frequencies for a Cen A and B are given 
in Table 12 

5. Comparison with previous studies 

We initially consider the calibrations made without seis- 
mic constraints and in particular the two recent de- 
tailed studies by Guenther & Demarque 1)2000(1 and Morel 
et al. pinnji . Contrary to our calibration, Guenther & 
Demarque (POOO ) did not use the masses of Pourbaix et 
al. (2002) which were of course not available at the time. 
However, they considered three different parallaxes which 
led to masses for the a Cen system very close to the ones 
of Pourbaix et al. I|2002|l . Indeed, they used masses con- 
tained between 1.0844 and 1.1238 Mq for a Cen A and 
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Fig. 8. Echelle diagram of a Cen A for the best model 
of the a Cen system (model M2), with a large spacing 
Al' ~ 105.5 /iHz. Open symbols refer to theoretical fre- 
quencies, while the filled circles correspond to the frequen- 
cies observed by Bouchy & Carrier (^2002 ,1 . Open circles are 
used for modes with 1 = 0, triangles for Z = 1, squares for 
I = 2 and pentagons for Z = 3 (see text for more details). 
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Fig. 9. Echelle diagram of a Cen B for the best model 
of the a Cen system (model M2), with a large spacing 
Av = 161.1 /zHz. Open symbols refer to theoretical fre- 
quencies, while the filled circles correspond to the frequen- 
cies observed by Carrier & Bourban 1)2003(1 . Open circles 
are used for modes with 1 = 0, triangles for Z = 1, squares 
for I = 2 and pentagons for I = 3 (see text for more de- 
tails). 

0.9017 and 0.9344 Mq for a Cen B (see Sect. [T] for more 
details). Given that these masses are similar to the ones 
we used for our calibration and that the input physics 
of both evolution codes is similar (mixing-length theory 
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Table 4. Low degree p-mode frequencies (in /iHz) for the best model (model M2) of a Ccn A and B. The observations 
are from Bouchy & Carrier l,2002j for a Cen A and from Carrier & Bourban 1,2003 1 for a Cen B. 







Observations 




Model M2 






n 


/ = 


I = 1 


I = 2 


I = 


I = 1 


I = 2 


I = 3 




15 






1833.1 


1717.8 


1764.7 


1814.0 


1857.8 




16 


1841.3 


1887.4 


1934.9 


1820.7 


1868.4 


1917.9 


1962.0 




17 




1991.7 


2041.5 


1924.2 


1972.4 


2022.9 


2067.7 




18 




2095.6 


2146.0 


2029.1 


2077.7 


2128.1 


2173.6 




19 
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2202.8 
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2279.1 


a Cen A 
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2339.4 


2385.2 




21 


2364.2 


2414.3 


2464.1 


2344.6 
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2444.8 


2491.2 
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2449.7 


2500.0 


2551.0 


2597.4 




23 


2573.1 


2625.6 




2555.5 


2605.8 


2657.3 


2704.3 




24 


2679.8 


2733.2 


2782.9 


2661.5 


2712.4 


2763.6 


2810.9 




25 


2786.2 


2837.6 


2887.7 


2767.5 


2818.6 


2870.3 


2917.5 




17 




3052.0 




2977.2 


3053.2 


3125.6 


3194.2 




18 








3136.7 


3213.2 


3286.2 


3355.1 




19 








3296.9 


3373.4 


3446.8 


3516.3 




20 




3532.9 


3607.2 


3457.2 


3534.0 


3607.5 


3677.4 




21 


3617.2 






3617.6 


3694.5 


3768.5 


3838.8 


a Cen B 


22 




3855.0 


3930.6 


3778.2 


3855.4 


3929.8 


4000.7 




23 


3937.6 


4016.3 




3939.2 


4016.8 


4091.5 


4163.0 




24 




4177.4 


4253.8 


4100.6 


4178.6 


4253.7 


4325.4 




25 








4262.6 


4340.6 


4415.9 


4488.2 




26 


4423.1 






4424.6 


4502.8 


4578.3 


4650.7 




27 


4585.0 






4586.7 


4665.0 


4740.8 


4813.5 



for convection and helium and heavy-element diffusion), 
one expects good agreement bctvirecn their results and 
ours. This is indeed the case. First of all, the initial he- 
lium aboundance of our models {Yi — 0.275 ± 0.010) is in 
good agreement with the value of about 0.28 obtained by 
Guenther & Demarque. Secondly, Guenther & Demarque 
found an age of 6.8 ± 0.8 Gyr if a Cen A does not have 
a convective core. As explained in Sect. 14.21 our model of 
a Cen A lies very close to the boundary between mod- 
els with and without a convective core, but does not ex- 
hibit a convective core. Thus, the age of 6.52 ± 0.30 Gyr 
of our M2 model is perfectly compatible with the results 
of Guenther & Demarque. Finally, Guenther & Demarque 
found that the mixing-length parameter of a Cen A is 
about 10% smaller that the one of the B component. 
However, they pointed out that this difference may not 
be significant, because the observables available at the 
time of their studies were not sufficient to strongly con- 
strain the models. Thanks to the asteroseismic results of 
Bouchy & Carrier (J002j and Carrier & Bourban (2003), 
we obtained a model of the a Cen system which is now 
firmly constrained. As a result, we confirm that a single 
mixing-length parameter cannot be found for both com- 
ponents of the system, as already suggested by many au- 
thors (Noels et al. lTMTl Fernandes & Neuforge UnHP and 
Guenther & Demarque 1^000,1 . Indeed, our results show 
that this parameter is about 8 % larger for a Cen B than 
for a Cen A, in perfect agreement with the value of about 
10% determined by Guenther & Demarque. Note that. 



at first sight, we would expect this parameter, if physi- 
cally meaningful, to be the same for all stars regardless of 
mass, composition and age. However, as shown by models 
with rotational mixing or mixing by magnetic instability, 
many more effects act, to the first order, as a change of the 
mixing-length parameter. Thus, differences in a = I /Hp 
must rather be considered as a measure of our ignorance 
in the internal stellar hydrodynamics, rather than only 
a difference stricto sensu of the mixing-length parame- 
ter. Note also that the mixing-length parameter of the A 
component is only slightly larger than the solar calibrated 
mixing-length parameter {aQ — 1.7). 

The calibration of Morel et al. (j2D00 ) was based on the 
masses obtained by Pourbaix et al. (jl999)l : Ma = 1.16 ± 
0.031 A/q and Mb = 0.97 ± 0.032 M©. These masses are 
larger than the ones considered by Guenther & Demarque 
(|2()()()ll and the new masses of Pourbaix et al. plH^ that 
we used for our models. Accordingly, the results obtained 
by Morel et al. H2()()()|l are quite different. For example, 
Morel et al. found an age of only 2.71 Gyr (for their model 
with the mixing-length theory and without overshooting). 
They also found values of the convection parameter that 
are almost equal for both stars. All these differences are 
mainly due to the different masses used. Note also that 
Morel et al. used the spectroscopic gravities to constrain 
the models instead of the luminosities derived from the 
photometry. This results in different error boxes in the 
HR diagram that have to be matched by the models and 
thus also inffuences the results of the calibration. 
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The first study of the a Cen system which took into 
account the asteroseismic measurements of the A compo- 
nent (Bouchy & Carrier OUn^ was made by Thevenin et 
al. I)2()02 V At first sight, one can think that our results 
are incompatible with this study, since Thevenin et al. 
pointed out that they were unable to produce a model of 
a Cen A and B compatible with the non-asteroscismic 
constraints and the data of Bouchy & Carrier (,20,02) us- 
ing the masses of Pourbaix et al. H2()()2|l . To understand 
why we are able to determine a model of the a Cen sys- 
tem which is in good agreement with all observables now 
available for these stars, one has to compare the non- 
asteroseismic constraints used by the two groups. As ex- 
plained in Sect. 12. 2l we adopted the same effective temper- 
ature as Thevenin et al. for a Cen B, while we increased 
the uncertainty on the effective temperature of a Cen A 
(5810 ± 50 K) in order to encompass the intervals of tem- 
perature used by Thevenin et al. and Thoul et al H2003p . 
In the same way, we adopted [Fe/H]A = 0.22 ± 0.05 dex 
for a Cen A and [Fe/H]B = 0.24 ± 0.05 dex for a Cen B; 
these values lie between the ones used by Thevenin et al. 
and those adopted by Thoul et al, with larger error boxes 
which appear to us more realistic in view of the different 
values found in the literature. Finally, Thevenin et al. de- 
termined the luminosities of both stars from the Geneva 
photometry using mean values with errors which do not 
take into account the quality of the night. We also used 
the Geneva photometric system, but calculated the mean 
magnitudes using the individual values listed in Table ^ 
and taking into account the quality of the night. As a 
result, the luminosities used by Thevenin et al. are very 
similar to our luminosities but with smaller uncertainties. 
It is important to underline these small differences, be- 
cause they can explain the different results found by the 
two studies. Indeed, our model of the a Cen system gives 
a luminosity of a Cen B which is in good agreement with 
our adopted observational constraints, but not with the 
small error box in the HR diagram adopted by Thevenin 
et al. Consequently, it is not surprising that Thevenin et 
al. did not find the solution we obtained. Moreover, the 
input physics of the evolution code used by Thevenin et al. 
differs from the input physics of our code. The main differ- 
ence concerns the modelisation of convection: Thevenin et 
al. used the Canute and Mazzitelli formulation while we 
used the mixing-length theory. These two different formu- 
lations can lead to significant different results, as shown 
by Morel et al. (|5000) . Another important difference con- 
cerns the hypothesis of a unique convection parameter for 
both stars. Indeed, Thevenin et al. assumed the convec- 
tion parameters to be the same for both components. As 
aforesaid, we found that a single mixing-length parame- 
ter is not applicable for both components of the system, 
in good agreement with Thevenin et al. 

We also notice that the model proposed by Thevenin 
et al. is not in good agreement with the new asteroseismic 
measurements of a Cen B (Carrier & Bourban 2003) due 
to theoretical large and small spacings that are higher than 
observed. Besides, Kervella et al. (|2()()3|l slightly changed 



the masse and the mixing-length parameter of the B com- 
ponent in order to obtain a model of a Cen B with a 
slightly larger radius than the one found by Thevenin et al. 
which is compatible with the interferometric radius. This 
larger radius for a Cen B certainly improves the agree- 
ment between their model and the observed mean large 
spacing. It is perhaps possible to obtain a model based on 
the masses derived by Thevenin et al. which is also com- 
patible with the observed mean small spacing by redoing 
the whole calibration in the same way as Thevenin et al., 
but with the asteroseismic measurements of a Cen B as 
additional constraints. However, we consider the observed 
masses as true observables which have to be reproduced 
by a consistent model of the a Cen system. 

Finally, we compare our results to the second calibra- 
tion of a Cen A and B which took into account the as- 
teroseismic observations of a Cen A (Thoul et al. 12003(1 . 
At first sight, one can think that the study by Thoul et 
al. is in good agreement with our work, since they de- 
termined a model of the a Cen system which reproduced 
the observed p-mode frequencies as well as the masses 
of Pourbaix et al. H2002|l . Moreover, they found an age 
of 6.41 Gyr which is in good agreement with our value of 
6.52 ± 0.30 Gyr. However, we have to notice that the non- 
asteroseismic observational constraints they used are quite 
different from ours. Indeed, Thoul et al. used the spec- 
troscopic surface gravities, while we used the luminosities 
derived from the photometry, bolometric correction and 
parallax. We chose to use the luminosities instead of the 
surface gravities to constrain the models, because the lu- 
minosities are determined with higher accuracy than the 
surface gravities. Consequently, our error boxes in the HR 
diagram are smaller than the ones adopted by Thoul et al. 
As a result, we see that our model of a Cen B lies within 
the error box adopted by Thoul et al., whereas the model 
of a Cen B determined by Thoul et al. is not included in 
our error boxes; this simply means that our model is also 
in perfect agreement with the surface gravities, whereas 
the solution obtained by Thoul et al. is not compatible 
with the observed luminosity of a Cen B. Another main 
difference between both studies concerns the input physics 
of the evolution code. Indeed, the evolution code used for 
our calibration includes a detailed treatment of the atomic 
diffusion of helium and heavy elements, contrary to Thoul 
et al. who neglect it. This fact can also explain certain dis- 
crepancies between the results found by Thoul et al. and 
our results, and in particular the fact that they obtained 
an identical mixing-length parameter for both stars. 

6. Conclusion 

The aim of this work was to determine the best model 
for the a Cen system using the Geneva evolution code 
including atomic diffusion. This model had to reproduce 
all observational constraints available for a Cen A and 
B, namely the masses, the luminosities, the effective tem- 
peratures, the metallicities, the radii and the low degree 
p-mode frequencies of both stars. First, we used all non- 
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asteroseismic constraints with one-sigma error boxes and 
found the solution t = 6.50 ± 0.20 Gyr, aA = 1.83 ± 0.10, 
aB = 1-99 ± 0.10, Yi ^ 0.275 ± 0.010 and {Z/X)i = 
0.0435 ± 0.0020. However, this model, which is in per- 
fect agreement with all non-asteroseismic measurements, 
was found to have large spacings Avq slightly larger than 
observed (105.9 ± 0.1 ^iRz instead of 105.5 ± 0.1 ^iHz for 
a Cen A and 161.7 ± 0.1 /xHz instead of 161.1 ± 0.1 ^iHz 
for a Cen B). These small discrepancies indicate that 
the radii observed by interferometry are slightly smaller 
than the ones deduced from asteroseismic measurements. 
Accordingly, we increased the error boxes on the ob- 
served radii from one to two sigma and found the fol- 
lowing solution: t = 6.52 ± 0.30 Gyr, ua = 1.83 ± 0.10, 
aB = 1-97 ± 0.10, Yi = 0.275 ± 0.010 and {Z/X)i = 
0.0434 ± 0.0020. Note that the radii of this solution are 
still in good agreement with the interferometric results of 
Kervella et al. PTMjl since they are only 0.4% and 0.8% 
larger than the observed radii for a Cen A and B respec- 
tively. 

To sum up, we point out that we obtained a model of 
the a Cen system which correctly reproduce all the nu- 
merous observational constraints now available for both 
stars. Thanks to the recent asteroseismic measurements 
of both components (Bouchy & Carrier I!j002l and Carrier 
& Bourban .2003 ) . this model is firmly constrained. We 
therefore conclude that asteroseismic measurements are 
needed to determine accurate stellar parameters of a given 
star as well as to test the physics used in the stellar evolu- 
tion codes. Indeed, non-asteroseismic constraints are not 
able to discriminate stellar models with similar positions 
in the HR diagram but different internal structures. It is 
however worthwhile to recall that asteroseismic measure- 
ments by themselves are not sufficient to obtain a reliable 
model of a given star nor to test the physics of stellar 
models. Indeed, the analysis of a Cen A and B shows that 
the combination of asteroseismic and non-asteroseismic 
constraints is the only way to obtain interesting and ac- 
curate results. Consequently, the determination of precise 
non-asteroseismic parameters like the luminosity, the ef- 
fective temperature or the metallicity of a given star is 
needed. For this purpose, the analysis of a visual binary 
system constitutes an ideal target for future asteroseismic 
measurements, since it is the guarantee to obtain reliable 
non-asteroseismic parameters and additional constraints 
as the same initial chemical composition and age for both 
stars. 
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